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This  paper  describes  a  highly  interactive  method  for  computer  visualization  of  the  basic  physical  process 
of  dispersion  and  mixing  of  fluid  elements  in  convection-diffusion  systems.  It  is  based  on  transforming 
the  vector  field  from  a  traditionally  Eulerian  reference  frame  into  a  Lagrangian  reference  frame.  Fluid 
elements  are  traced  through  the  vector  field  for  the  mean  path  as  well  as  the  statistical  dispersion  of  the 
fluid  elements  about  the  mean  position  by  using  added  scalar  information  about  the  root  mean  square 
value  of  the  vector  field  and  its  Lagrangian  time  scale.  In  this  way,  clouds  of  fluid  elements  are  traced 
not  just  mean  paths.  We  have  used  this  method  to  visualize  the  simulation  of  an  industrial  incinerator 
to  help  identify  mechanisms  for  poor  mixing. 
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1  Introduction 


Many  computational  simulations  in  engineering  involve  solutions  of  simultaneous  vector  and  scalar  fields 
for  convection  and  diffusion  equations  in  three  dimensional  space.  The  physical  process  of  dispersion 
and  mixing  within  such  systems  is  often  what  the  engineer  is  trying  to  understand  when  performing 
the  numerical  simulation.  Current  computer  visualizations  of  calculated  scalar  and  vector  fields  do  not 
provide  sufficient  insight  to  immediately  conceptualize  the  underlying  dispersion  and  mixing  process 
inherent  in  these  systems.  For  example,  in  computational  fluid  dynamics  problems  involving  the  mixing 
and  reaction  of  three  input  streams  (A,  B  and  C),  the  user  would  like  to  know  how  quickly  stream  A 
disperses,  where  it  is  going,  and  if  and  where  it  mixes  with  streams  B  and  C  which  are  also  convecting 
and  dispersing.  Schemes  to  visualize  the  scalar  concentration  of  stream  A  tell  something  about  the 
extent  of  mixing  but  nothing  about  where  A  came  from  or  with  which  stream  (B  or  C)  it  has  mixed. 
Particle  tracing  methods  can  give  an  idea  of  where  the  mean  flow  path  is  found  for  a  representative 
massless  fluid  element  starting  at  any  point  (i.e.  in  stream  A)  but  nothing  about  its  dispersion  rate  or 
its  interaction  (mixing)  with  other  streams. 

In  this  paper  a  highly  interactive  visualization  scheme  is  presented  which  attempts  to  allow  a  compu¬ 
tational  fluid  dynamics  user  to  visualize  the  basic  physical  process  of  dispersion  and  mixing  rather  than 
just  the  vector  or  scalar  values  computed  by  the  simulation.  The  motivation  for  this  approach  is  the 
hope  that  if  the  user  can  gain  a  greater  insight  into  convection  and  diffusion  through  the  graphical  visu¬ 
alization,  then  the  original  simulation  will  be  of  greater  utility  in  solving  practical  engineering  problems. 
As  shown  in  our  test  results,  this  visualization  method  allows  us  to  identify  faults  in  the  design  of  of  an 
industrial  incinerator. 

The  method  is  based  on  transforming  the  vector  field  from  a  traditionally  Eulerian  (fixed)  reference 
frame  into  a  Lagrangian  (moving)  reference  frame.  Mean  fluid  elements  are  traced  through  the  computed 
vector  field  (as  in  traditional  particle  tracing  methods).  The  statistical  dispersion  of  the  fluid  elements 
about  the  mean  position  is  also  computed  by  using  added  scalar  information  about  the  root  mean  square 
value  for  the  vector  field  and  its  Lagrangian  time  scale.  These  scalar  values  can  usually  be  computed 
easily  from  existing  information  in  computational  fluid  dynamics  codes  (i.e.  turbulence  models,  etc.). 
In  this  way  clouds  or  plumes  of  fluid  elements  are  traced  not  Just  mean  paths. 

The  cloud  surface  is  represented  as  a  triangular  mesh  which  can  be  efficiently  rendered  by  conventional 
graphics  hardware.  The  surface  is  colored  by  using  any  scalar  field  such  as  temperature,  pressure,  etc. 
selected  by  the  user.  During  a  visualization  session,  the  user  can  interactively  pick  multiple  three- 
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dimensional  positions  in  the  data  domain  and  see  how  flows  starting  from  these  positions  develop  and 
mix  together  in  some  way.  Local  twisting  or  deformation  can  be  visualized  by  displaying  ribbons  along 
the  mean  trajectory.  The  interactive  cloud  tracing  scheme  presented  in  this  paper  complements  existing 
visualization  techniques. 

2  Flow  Visualization 

As  both  the  cost  of  experimental  flow  visualization  and  the  scale  of  the  problems  that  scientists  can  sim¬ 
ulate  continue  to  increase,  flow  visualization  using  computer  graphics  techniques  becomes  an  essential 
and  popular  mechanism  for  scientific  researchers  to  understand  the  results  of  fluid  dynamics  simulations 
in  fields  such  as  aeronautics,  transportation  design,  weather  forecasting,  oceanography,  combustion  en¬ 
gineering,  etc.  Well  known  flow  visualization  methods  making  use  of  computer  graphics  include  arrow 
fields,  contouring,  streamlines,  particle  tracing,  etc.  To  date,  many  elaborate  methods  for  graphical 
flow  visualization  have  been  developed,  but  none  of  them  has  emerged  as  being  the  standard  due  to  the 
complexity  of  three-dimensional  flow  fields. 

Stream  surfaces/ribbons  [9]  provide  information  about  twisting  along  particle  paths.  The  stream 
polygon  method  [14]  displays  local  deformation  due  to  normal  and  shear  strain  and  rigid  body  rotation. 
Silver  et  at.  [15,  18]  describes  how  to  visualize  vortex  structures  by  employing  so  called  visiometric 
techniques  including  thresholding,  object  isolation,  abstraction  and  data  juxtaposition. 

The  concept  of  using  many  moving  particles  to  create  the  illusion  of  fluid  movement  originated  from 
Reeves  [13]  Spot  noise  [17]  is  the  use  of  stochastic  texture  (overlapping  particles)  for  the  visualization  of 
scalar  and  vector  fields  over  surfaces.  A  spot  is  a  two-dimensional  pattern  like  a  disk,  a  square  or  pair 
of  crossed  lines.  The  attributes  of  spots  such  as  size,  shape,  direction,  and  edge  are  varied  according 
to  local  scalar  and  vector  flow  fields  resulting  in  directional  information,  which  provides  a  qualitative 
visualization  of  the  data.  Particle-based  rendering  is  usually  computationally  expensive.  Crawfis  and 
Max  [7,  1 1]  develop  an  efficient  technique  integrating  a  fine-grmned  texture  representation  of  vector  flow 
field  into  a  direct  volume  rendering  method  to  visualize  vector  and  scalar  fields  simultaneously.  Gelder 
and  Wilhelms  [16]  present  a  highly  interactive  flow  animation  technique  also  based  on  the  concept  of 
particle  systems.  Particles  are  placed  randomly  in  the  flow  field.  Each  particle  path  is  represented  as  a 
polyline.  Then  particle  motion  is  created  by  using  a  color-table  animation  technique. 

Critical  points  are  locations  in  the  flow  field  where  the  velocity  field  vanishes.  These  points  can  be 
characterized  according  to  the  nearby  integral  curves.  Critical  points  in  a  flow  field  along  with  the 


2 


integral  curves  connecting  them  may  completely  describe  the  topology  of  the  vector  field.  Much  work 
has  been  done  in  the  visualization  of  two-dimensional  flow  topology  and  three-dimensional  separation 
flows  near  the  surface  of  a  body  [4,  6,  8].  The  integration  of  critical-point  analysis  into  a  graphical  flow 
visualization  process  also  provides  a  powerful  divide-and-conquer  approach  to  data  exploration  since  the 
topological  information  can  be  used  to  partition  a  large  and  complex  flow  field  into  well-defined  regions. 

Virtual  smoke  is  an  interactive  flow  tracing  technique  that  we  developed  previously  [10].  It  simulates 
the  use  of  colored  smoke  for  experimental  gaseous  fluid  flow  visualization.  It  is  based  on  direct  volume 
rendering  so  realistic  three-dimensional  visualization  of  the  flow  can  be  produced.  In  particular,  it  can 
animate  the  dynamic  behavior  of  steady-state  or  instantaneous  flow  fields.  While  the  virtual  smoke 
method  can  display  realistic  flow  patterns  in  motion,  its  diffusive  component  is  less  accurate  and  it  in 
fact  can  produce  diffusive  flow  artificially. 

The  cloud  tracing  method  presented  in  this  paper  is  a  revision  of  virtual  smoke  to  make  visualizations 
that  are  more  physically  accurate.  The  next  section  gives  a  brief  description  of  the  physical  model  and 
the  related  mathematical  equations  behind  the  cloud  tracing  visualization  method.  Then  the  graphics 
algorithm  used  to  realize  the  cloud  tracing  visualization  is  described  in  Section  4.  We  have  applied 
this  tracing  technique  to  a  data  set  from  the  simulation  of  an  industrial  incinerator  and  the  results  are 
reported  in  Section  5. 
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3  A  Dispersion  Model 

The  dispersion  of  particles  or  fluid  elements  by  turbulent  motion  has  been  an  active  field  of  research  in 
recent  years.  Practical  interests  are  widespread  and  include  the  combustion  of  particulate  and  droplet 
fuels,  the  pneumatic  transport  of  particles,  the  mixing  and  settling  of  industry  emissions,  particle  de¬ 
position,  electrostatic  precipitation,  and  particle  fluid  separations.  Practical  combustion  processes  like 
spray  combustion,  pulverized  coal  and  coal  slurry  combustion,  fluidized  beds,  sorbent  injection,  and 
hazardous  waste  incineration,  introduce  the  fuel  as  particles,  liquid  droplets,  or  slurries  into  turbulent 
environments. 

.Mixing  and  dispersion  of  gaseous  streams  involves  fluid  elements  that  disperse  like  particles  but  without 
any  added  inertial  effects  of  particles  with  drag.  The  equations  used  here  were  developed  by  Batchelor 
[Ij.  They  have  been  extended  and  used  for  turbulent  dispersion  of  inertial  particles  in  a  model  for 
the  stochastic  transport  of  particles,  referred  to  as  the  STP  model  [2].  This  approach  is  a  stochastic 
Lagrangian  scheme  that  can  accurately  and  efficiently  compute  fluid  dispersion  in  complex  turbulent 
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Figurt*  1 :  A  La^rangian  View  of  Dispersion  in  a  One-dimensional  Flow 


flows. 


In  this  approach,  the  time  evolution  of  a  probability  density  function  for  fluid  particle  position  is 
modeled.  Figure  1  illustrates  pictorially  the  development  of  a  probability  density  function  for  particle 
position  for  an  initially  concentrated  ensemble  of  fluid  particles  released  from  a  source  near  the  origin 
in  a  two  dimensional  coordinate  system.  The  circles  circumscribe  the  region  of  the  flowfield  included 
within  one  standard  deviation  of  the  mean  fluid  particle  position,  which  is  at  the  center  of  each  circle. 
A  brief  mathematical  description  of  the  approach  is  given  in  the  following  section. 

The  path  of  an  individual  fluid  element  is  given  by 


x(t)  =  [  V{T)dT  -I-  x(0) 
Jo 


:i; 


where  V  is  fluid  velocity  and  x  is  particle  position. 

Then,  for  a  stationary  velocity  field,  the  variance  of  the  probability  density  function  for  fluid  particle 
position  can  be  written  as 


=  2  f  (t-  r)«(r)dr 
Jo 


(2) 
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R(t)  is  the  fluid  particle  autocorrelation  function  and  takes  the  form 

R{r)  =  (3) 

where  Vms  is  the  root  mean  square  fluid  velocity  fluctuation,  and  Ti  is  the  Lagrangian  integral  time 
scale  which  is  the  average  persistence  of  the  turbulence  activity  at  a  point.  The  root  mean  square  gas 
velocity  fluctuations  are  assumed  to  be  isotropic  and  are  calculated  from  the  turbulent  kinetic  energy, 
k,  as 

0  2k 

Vrrr..  =y 

Tl  comes  from  the  description  of  the  turbulent  field.  For  example,  if  a  k-e  model  is  used, 

^  _  Cjk^^ 

I L  —  1  • 

where  is  a  constant  and  £  is  the  computed  dissipation  rate  of  turbulence. 

4  Visualization  of  Cloud  Tracing 

I'nlike  traditional  computer  graphics,  the  goal  of  visualization  is  insight,  not  reabstically  rendered  images. 
The  focus  of  the  following  discussion  is  on  a  surface  rendering  technique,  which  allows  the  user  to  depict 
how  fluid  elements  retain  certain  properties  in  a  highly  interactive  manner.  The  surface  normally  appears 
as  a  growing  cone  which  is  deformed  due  to  fluctuating  local  vector  fields.  By  using  surface  rendering,  we 
sacrifice  realism  and  hope  that  higher  interactivity  and  flexibiUty  can  better  assist  the  user  in  obtaining 
insight  into  the  flow  fields.  As  described  later,  more  reaUstic  visualization  of  the  flow  can  still  be  obtained 
by  using  a  direct  volume  rendering  technique  similar  to  virtual  smoke.  By  any  means,  the  fluid  particle 
dispersion  is  calculated  statistically  and  the  diffusive  behavior  of  the  flow  is  described  more  accurately. 

4.1  TVacing 

The  cloud  tracing  starts  from  a  three-dimensional  position  P(x,y,z)  specified  by  the  user  using  a  three- 
dimensional  cross  hair.  At  this  position,  V ,  Vrm*  and  Tl  are  trilinearly  interpolated  from  the  neighboring 
eight  voxels.  As  in  traditional  fluid  particle  tracing,  knowing  V ,  P  and  a  user-specified  time  increment  dt, 
the  mean  particle  trajectory  T  can  be  calculated  using  a  numerical  integration  method  like  the  two-step 
Runge-Kiitta  time-advance  method  such  that: 

P*  =  P,+dtV{P,) 

P,+,  =  P-  +^-dt{V{P-)-V{P,)) 
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Figure  2:  Dispersion  Area  as  a  Function  of  Variance  within  3  Standard  deviation. 

Thus  T  contains  a  sequence  of  points  Pi{x,  y,  z)  along  the  trajectory  calculated  where  Po  =  P.  At  each 
mean  point  Pj ,  we  then  calculate  the  dispersion  of  the  fluid  elements  about  P,  by  using  scalar  quantities 
Vrma  &nd  Tj,.  Equation  2  is  used  for  this  calculation.  ff{t)  defines  the  spatial  extent  of  the  probability 
density  function  for  particle  position  at  successive  and  equally  spaced  residence  times.  Intuitively,  a{t) 
defines  a  circle  circumscribing  the  region  of  the  flow  field  included  within  one  standard  deviation  of  the 
mean  particle  position.  That  is,  Pj  is  the  center  and  a{t)  is  the  radius  of  the  circle.  Figure  2  shows  the 
dispersed  area  as  a  bivariate  normal  distribution  function  of  the  variance  within  three  standard  deviation 
of  the  mean  particle  position.  The  integration  is  done  by  polynomial  interpolation.  Using  the  TVapezoid 
method.  Equation  2  becomes 

<r(t)*  =  2  f  (t-  T)R{T)dr 
Jo 

=  ^{tR(0)  +  2j2it-ih)Riih)) 

f=l 

where  h  =  ^  and  n  is  the  number  of  points. 

In  actual  implementation,  the  trajectory  is  calculated  incrementally  along  with  the  calculation  of 
the  variance  of  the  probability  density  function.  For  each  pair  of  consecutive  points  Pi  and  P.+i,  a 
triangle  mesh  connects  the  circumference  of  the  two  circles  defined  at  the  corresponding  points.  The 
orientation  of  a  circle  is  defined  by  the  velocity  at  the  corresponding  point.  Applying  a  triangle  mesh 
along  the  trajectory  in  this  fashion  generates  a  cone-shape  surface  which  represents  the  dispersion  along 
the  trajectory. 


4.2  Surface  Construction 

The  way  the  surface  is  constructed  is  similar  to  the  traditional  method  for  extracting  surfaces  between 
consecutive  contours  [12].  However,  there  is  no  ambiguity  in  our  case  since  the  consecutive  slices  to 
be  connected  are  always  perfect  circular  planes  with  no  holes.  Our  algorithm  uses  a  standard  circular 
plane  which  is  transformed  from  point  to  point  according  to  the  local  velocity  V  =  {u,v,w).  we  define 
the  standard  circular  plane  on  the  xy  plane  and  centered  at  the  origin  of  the  corresponding  Cartesian 
coordinate  system.  In  addition,  the  derivation  of  a  transformation  matrix  for  the  rotation  of  a  vector 
about  a  general  axis  through  the  origin  that  is  described  in  Faux’s  book  [5]  can  be  used  for  our  purpose. 
Because  the  orientation  of  the  standard  plane  is  (0, 0, 1 .0),  the  transformation  matrix  is  greatly  simplified 
as  follows: 

v^{\  —  w)  +  w  —uv(l  —  w)  sin(ff)u  0 
—  ut;(l  —  tn)  u^(l— 2)  +  iu  sin{0)v  0 
-sin{$)u  -sin(ff)v  w  0 

X  y  z  1 

where  «,  v,  and  w  are  normalized,  0  =  y/u^  +  which  is  the  angle  between  (0,0, 1.0)  and  (u,  v,  w),  and 
X,  y,  and  z  are  the  coordinates  of  P,.  Figure  3  shows  the  standard  plane  and  a  sequence  of  transformations. 
Note  that  only  the  points  C,(n)  approximating  the  circle  are  transformed.  The  number  of  points  used 
is  selected  by  the  user.  Alternatively,  the  number  of  points  can  grow  adaptively  according  to  the  radius 
of  the  circle;  however,  care  must  be  taken  to  ensure  continuity  of  the  triangle  mesh.  Note  that  a  4-point 
model  would  resemble  the  stream  polygon  visualization  [14]  but  is  not  sheared.  The  surface  normal  is 
approximated  by  calculating  C,(n)  -  P,  at  each  C,(n).  The  surface  is  hardware  Gouraud  shaded.  The 
interpolated  scalar  value  at  each  triangle  point  is  mapped  to  the  color  selected  by  the  user.  Figure  4 
shows  the  wireframe  representation  of  two  cloud  surfaces. 

In  addition,  a  ribbon  tracer  can  be  made  by  consistently  selecting  two  diagonally  opposite  points  on 
each  circle  and  then  constructing  a  surface  using  those  points  (spinwise  edges)  as  shown  in  Figure  5.  A 
cross  ribbon  tracer  can  be  made  by  selecting  two  sets  of  such  diagonal  points  on  each  circle.  However, 
the  two  lines  crossing  the  two  pairs  of  points  on  each  circle  should  be  orthogonal.  Besides  being  used  to 
examine  local  twisting  and  deformation,  ribbons  may  also  help  observe  flow  mixing.  The  cloud  tracing 
technique  described  above  has  been  implemented  in  Motif  and  GLX  for  running  on  Silicon  Graphics 
workstations.  Motif  is  used  to  implement  the  user  interface.  GLX  which  allows  GL  and  X  mixed¬ 
mode  programming  is  used  to  implement  3D  graphics.  On  a  workstation  supporting  alpha  blending, 
transparent  surfaces  can  be  used  to  derive  a  more  intuitive  visualization. 
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v=(  u,  V,  w  ) 

Figure  3:  Standard  Plane  and  Transformations. 


A  realistic  volume  rendered  image  of  cloud  tracing  can  still  be  obtained  after  the  traces  are  made 
by  using  the  technique  described  in  [3,  10].  Essentially,  the  mean  particle  paths  P,{t)  are  treated  as 
seedlings  and  (T(t)  is  used  as  the  radius  to  the  corresponding  seed.  The  result  is  an  improved  virtual 
smoke  which  accurately  visualizes  the  convection  as  well  as  diffusive  flow.  Nevertheless,  interactive  direct 
volume  rendering  is  still  limited  to  a  fixed  view  position. 

5  Test  Results 

A  computer  simulation  of  an  organic  waste  incinerator  was  performed.  In  this  cylindrical  furnace  the 
waste  is  injected  into  the  combustion  chamber  through  28  burners  located  around  the  perimeter  and  near 
the  base  of  the  furnace.  Each  burner  has  two  concentric  annuli;  the  waste  stream  is  injected  through 
the  center  and  air  through  the  outer  annulus.  The  objective  of  the  incinerator  design  is  to  accomplish 
complete  mixing  and  reaction  of  the  waste  and  air  streams  as  rapidly  as  possible.  The  objective  of  the 
visualization  is  to  explore  the  success  or  failure  of  the  design  for  accomplishing  this  rapid  mixing.  In  this 
test,  because  of  the  symmetry  of  the  furnace  design  and  the  area  of  interest,  only  one  side  of  the  lower 
half  of  the  furnace  is  visualized.  This  portion  of  the  furnace  contains  100,800  data  points  in  cylindrical 
coordinates. 


Plate  1  shows  traditional  fluid  particle  traces  originating  in:  1)  the  center  of  the  waste  stream  and  2) 
the  renter  of  the  air  stream,  for  one  of  the  28  burners.  Although  the  fluid  path  is  shown  in  this  method,  it 


Figurp  4:  Cloud  Constnirtion. 


Figure  5:  Ribbon  Construction. 


gives  little  information  on  the  extent  of  mixing  between  the  two  streams.  If  we  track  more  fluid  elements 
we  still  do  not  gain  any  better  insight  into  the  path  of  the  cloud  since  the  paths  that  are  traced  only 
represent  the  mean  location  for  any  group  of  fluid  elements  released  from  the  initial  location.  The  actual 
mixing  includes  the  statistical  dispersion  (convection  and  diffusion)  of  the  entire  cloud  released  from  the 
specified  location. 

Plate  2  shows  the  cloud  tracing  performed  by  the  technique  discussed  in  this  paper.  The  surface 
around  each  mean  path  encloses  one  standard  deviation  (68.27%)  of  each  cloud.  Now  the  extent  of 
mixing  of  these  two  clouds  is  apparent.  Plate  3  shows  a  view  that  is  rotated  to  look  from  the  top  of 
the  incinerator.  Because  this  furnace  was  designed  with  so  many  burners  and  with  both  air  and  waste 
in  each  burner,  it  was  expected  that  mixing  would  be  accomplished  rapidly.  Both  experimental  testing 
and  numerical  simulation  showed  otherwise.  The  objective  of  the  simulation  was  to  see  why  mixing  was 
not  as  rapid  as  was  expected.  Simple  color  shading  of  the  waste  concentration  at  the  exit  plane  of  the 
furnace  showed  high  waste  concentrations  near  the  center  of  the  furnace  and  high  concentration  of  air 
at  the  outer  perimeter  of  the  furnace.  The  cloud  tracing  visualization  identified  the  reasons  for  this. 
The  two  clouds  shown  in  Plate  2  and  3  show  that  the  waste  stream  from  each  burner  penetrates  further 
into  the  center  of  the  cylindrical  furnace  before  turning  upward  by  the  mean  flow  from  below;  while  the 
air  stream  on  the  outside  of  each  burner  turns  upward  much  sooner  and  stays  at  the  outer  perimeter  of 


the  furnace.  The  cloud  surface  is  shaded  by  mapping  O2  concentration  to  colors;  low  concentration  to 
white-blue  and  high  to  red-purple. 

The  rendering  of  the  dispersing  cloud  as  a  colored  surface  makes  real  time  interactivity  possible. 
Interactivity  of  the  visualization  scheme  is  necessary  for  the  analysis  of  complex  mixing  chambers.  The 
spatial  relation  of  the  dispersing  stream  tube  can  only  be  appreciated  when  the  rendered  image  can  be 
rotated  in  real  time.  The  scalar  field  is  mapped  to  color  on  the  tube  surface.  Real  time  rotations  allow 
a  full  exposure  of  the  colored  surface  to  the  observers’  view.  The  color  of  the  surfaces  in  Plate  2  and  3 
is  representative  of  the  air  concentration.  Thus  the  visualization  indicates  simultaneously  that  the  two 
jets  shown  do  not  mix  readily;  but  instead,  air  from  other  surrounding  jets  tends  to  stay  at  the  outer 
radius  of  the  incinerator.  Plate  4  shows  a  surface  enclosing  2.0  standard  deviations  (95%)  for  each  of 
the  clouds  traced.  Only  the  outer  portions  of  these  two  representative  clouds  begin  to  mix  by  the  end 
of  this  furnace. 

This  visualization  scheme  has  been  able  to  uncover  mechanistic  reasons  for  the  poor  mixing  in  an 
industrial  incinerator.  The  poor  mixing  was  not  expected  from  an  engineering  analysis  of  the  initial 
design.  The  reasons  for  the  poor  mixing  were  not  uncovered  from  experiments  performed  with  the 
incinerator  or  from  the  numerical  simulation  of  the  furnace  without  this  visualization  tool. 

6  Conclusions 

Interactivity  is  the  key  to  successful  scientific  visualization.  We  have  introduced  a  highly  interactive 
cloud  tracing  technique  for  the  visualization  of  simultaneous  vector  and  scalar  fields  in  convection- 
diffusion  systems.  A  dispersion  model  for  fluid  elements  is  integrated  into  the  visualization  process  to 
more  accurately  portray  the  physical  process  of  dispersion  and  mixing  within  such  system.  The  use  of 
a  surface  rendering  technique  allows  the  user  to  move  the  cloud  tracer  in  real  time.  By  using  a  data  set 
from  the  simulations  of  an  industrial  incinerator,  we  have  demonstrated  that  this  cloud  tracing  technique 
is  an  effective  analysis  tool  for  practical  engineering  problems.  It  gives  visualization  to  the  mechanism 
of  gaseous  mixing  not  just  the  result  of  mixing. 

Future  work  includes  a  better  user  interface  for  locating  seed  points,  the  integration  of  volume  ren¬ 
dering  for  viewing  more  than  the  cloud  surfaces,  tracing  of  particle  cloud  and  enhancement  for  handling 
unsteady  flows.  For  unsteady  flows,  clouds  must  evolve  simultaneously  and  not  sequentially  as  in  our 
current  implementation. 
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